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Abstract 

The quantum mechanical description of the world leads to predictions that defy our daily experi- 
ence when dealing with macroscopic objects. Schrodinger's famous cat illustrates most graphically 
the paradoxical situations that one encounters when extrapolating the quantum rules to systems 
possessing macroscopically distinguishable states. Leggett and Garg formulated this question qual- 

O" 1 



itatively by deriving a series of inequalities based upon the premises of macroscopic realism 



qual- 



which have been recently tested within several microscopic systems |3H7( and micrometer-sized 
t^. . superconducting circuits [8j. Here we present an experimental violation of Legget-Garg inequali- 



ty | ties (LGI) with millimeter-sized crystals whose state can be thought of as a cat of cats. Through 

OV 

o 



■ the reversible transfer of polarization state of single photons and the collective excitation states of 



crystals, we prepared superposition states of two macroscopic crystals separated by a distance in 
the mm range. The phase relation between the states in the superposition is governed by the dy- 
namics of a large number of individual atoms and can be controlled with a frequency detuned and 
polarization dependent atomic frequency comb 9l4l3l|. The quantum character of the considered 



global state is therefore governed by a macroscopic quantity, which allows for an additional level 
of complexity in the system. The observed violation of a LGI confirms the persistence of quantum 
coherence effects within macroscopic objects which can be seen with a naked eye. 
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Ever since the birth of quantum mechanics (QM) it has been stressed the difficulty to 
reconcile the behavior of quantum particles and our intuitive experience when dealing with 
macroscopic objects, which should occupy definite states at all times and independently of 
the observers. Some physicists, with Niels Bohr as a precursor, argued that the solution 
to this problem requires dividing the world into a microscopic domain described by the 
quantum formalism and the macroscopic level of measuring apparatus, which is postulated 
to behave classically. However, this division unavoidably carries a degree of arbitrariness, as 
illustrated by the famous Schrodinger cat gedanken experiment [yj], where the superposition 
state of an atom is shown to be transferable, in principle, to a superposition state of two 
macrocopically distinct states. 

Following the argument presented in the 1935 Einstein-Podolsky-Rosen seminal paper 



15] , Bell derived an experimentally testable inequality, in which the premises of local realism 



yield statistical predictions that disagree with those of quantum mechanics 16|. Numerous 
experiments have confirmed the validity of the QM predictions in spatially separated mi- 
croscopic systems, with various loop-holes regarding experimental deficiencies having been 
closed individually Ul9|. A loophole-free experimental violation of Bell inequality would 
lead to the sacrifice of either the locality or the realism assumption. On the other hand, 
in contexts where locality is not relevant, the focus is on the realism assumption at a level 
much closer to our everyday perception. Unlike the Bell inequality, which involves statisti- 
cal correlations between states of two spatially separated systems, LGI are concerned with 
the correlations of the state of a single system at different times (autocorrelations) 1]. In- 
troduced as the Bell inequalities in time, the violation of LGI excludes a hidden-variable 
description based upon the assumptions of macroscopic realism (a macroscopic system must 
at any time be in a definite one of its macroscopically distinct states) and the possibility 
to perform noninvasive measurements (measurements do not influence the actual state or 
the subsequent system dynamics of the system). LGI therefore provide a criterion to char- 
acterize the boundary between the quantum and the classical domains and could witness 
the persistence of coherence effects beyond the closed systems paradigm 20]. Testing the 
validity of standard quantum theory over larger and larger scales and increasing levels of 
complexity is important for assessing the relevance of genuinely quantum traits beyond mi- 
croscopic scales [21)] as well as analyzing the feasibility of unconventional decoherence models 



like those invoking spontaneous wave function collapse 22N25]. 
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Until now, the experimental violation of LGI has been performed with microscopic ob- 
jects, such as photons (3), Q], electron or nuclear spins js-7| and the micrometer-sized su- 
perconducting 'transmon' system In practice, the main experimental challenge comes 
from the implementation of truly non-invasive measurements. However, the use of weak 
measurements 8] and schemes involving ancillary systems [7| have recently allowed for the 
first experimental results in the microscopic domain. To avoid the need of performing mea- 
surements at intermediate times, as required in the original formulation of the inequalities, 
it has been suggested that the condition of noninvasive measurements can be replaced by 



the weaker assumption of stationary correlations |26|, [27|]. The resulting formulation is then 
easily testable by using projective measurements. Within this framework, inequalities are 
derived from the premises of: 

1) Macroscopic realism: A macroscopic object which has two available macroscopically dis- 
tinct states is at any time t in one of these states, i.e., the dichotomic observable M(t) 
always produces the value +1 or -1. 

2) Stationary: The conditional probability Qij(t\,t%) to find a system in state j at time t%, 
if it is in state % at time t\ only depends on the time difference {ti — t\). 



LG-like inequalities can then be formulated as follows 
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231, 



K = K(0, 2t) - 2K(0, t) > -1, 



K 



(1) 



K(0,2t) + 2K(0,t) > -1, (2) 

where K(ti,t2) denotes the autocorrelation function defined as K{t\,tz) = (M(ti)M(i 2 ))- 
This inequalities can be shown to be violated by quantum mechanical unitary dynamics. It 
should be noted that, in principle, the assumption of stationarity narrows down the class 
of macrorealist theories which are put to the test. However the stationarity assumption 
can lead to testable inequalities, on the one hand, and on the other hand, we should stress 
that it is an assumption that could be tested independently, as opposed for instance to the 
situation encountered with standard Bell inequalities, where the so called "fair sampling 

]. As a result, if the stationarity condition does hold for the 



assumption" is not testable 
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consider experimental set up, we argue that the inequalities of the form Eqs.f lTf2"j) provide a 
fair test for a broad class of macrorrealistic models in that set up. 

Here we present an experimental violation of this form of a LGI with macroscopic crystals. 
The validity of the stationary assumption in this system is experimentally verified with an 



(a) (b) (c) 




FIG. 1: (a). The H (blue) and V (red) parts of the atomic frequency comb (AFC) are prepared 
with a periodicity of A and a relative detuning of 5. (b). The pulse sequence for AFC preparation 
used in the experiment. The frequency of pump light are scanned in each cycle while the intensities 
of H (blue solid line) and V (red dashed line) polarized light are modulated periodically, (c). The 
readout efficiency at different times. 

independent test. The physical system undergoing evolution is composed of about 10 10 
ions, which are spatially distributed into two, 2.4 mm separated crystals. The experimental 
sample, which was presented in [ll] to demonstrate a memory hardware, is composed of two 
pieces of Nd 3+ :YVC>4 crystals (doping level 10 ppm) sandwiching a half- wave plate (HWP). 
The sizes of the two crystals are nearly equal with 1.4±0.01 mm length along the a-axis. 
The 4 ig/2 — > 4 F 3 /2 transition of Nd 3+ at about 879.705 nm in the Nd 3+ :YVC>4 crystal shows 



strong absorption of H polarized photons and little absorption of V polarized photons [29 ]. 
Here, H(V) denotes horizontal (vertical), which is denned to be parallel (perpendicular) to 
the crystal's c-axis. By using a 45° HWP, sandwiched between two parallel Nd 3+ :YVC>4 
crystals, the sample shows nearly equal absorption depth for H and V polarized photons. 
We utilize the atomic frequency comb (AFC) technique to realize the reversible transfer 



of photonic polarization state and atomic state j9l— Il3j] . By independently preparing two 
frequency detuned AFC in the two crystals, we obtained a polarization dependent AFC, 
which enable us to initialize, control and analyze the collective excitation states of the 
crystals. 

The AFC protocol requires a pumping procedure to tailor the absorption profile of an 
inhomogeneously broadened solid state atomic medium with a series of periodic and narrow 
absorbing peaks separated by A (see Fig. 1). The single input photon is then collectively 
absorbed by all the atoms in the comb. The atomic state with and without the photon 
excitation can be represented by |e)jv = Yu7 Cje~ tkZi e l27r5:it \gi gAr) |9| and \g)N = 
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\gi ■ ■ ■ gj ■ ■ ■ g N ), respectively. Here N is the total number of atoms in the comb; \gj) and \ej) 
represent the ground and excited states, respectively, of atom j\ Zj is the position of atom j; 
k is the wavenumber of the input field; Sj is the detuning of the atom with respect to the laser 
frequency and the amplitudes Cj depend on the frequency and on the position of atom j. Due 
to the periodical structure of AFC, Sj ~ rrijA with rrij of integer. With the two AFC of the 
two crystals prepared with the same periodicity A, when the input photon's polarization is 
chosen as H + e tip0 V, the atomic state of the two crystals will be excited simultaneously, and 
can be represented by \&(t) = l/^/2(\a) + \b)e l ^ 27TSt+v>0 ^) where \a) and \b) are the excitation 
states of atoms in the first and second crystals, respectively. 5 is the frequency detuning 
between the two AFC, as shown in Fig. 1(a). The two crystals are actually in an entangled 
state which can be rewritten as l/v / 2(|e)Ari|g)AT2 + |<?)Ari|e)Ar2e 4 ( 27n5 * +¥ ' )). This entanglement 
has been already demonstrated in a non-polarizing configuration 30] . 

Now we define the observable M(t) as the superposition state of the collective excita- 
tion in the two crystals, where the linear state \D) = l/v2(|a) + \b)) is used as |0) with 
eigenvalue equal to +1 and the orthogonal state \A) = l/v2(|a) — \b)) is used as |1) with 
eigenvalue equal to -1. In which sense can the states \D) and \A) be branded as 'macro- 
scopically distinct'? Despite both states involve delocalized excitations over a mm scale, 



attending to the original formulation of Leggett [2j, |3jj , the considered states do not have 
a sufficiently high degree of dis connectivity and therefore their macroscopicity would be ar- 
guable. However, and that is a crucial part of our proposal, the two states \D) and \ A) differ 
by a collective phase shared by ~ 10 10 ions and it is this phase which dictates whether or 
not the system behaves as a quantum coherent superposition or as a (classical) statistical 
mixture. We should emphasize that unlike the circulation quantization in superfluid helium 
or Josephson effect in superconductors, the acquired phase difference between states \a) and 
| b) cannot be simply interpreted as an accumulated microscopic quantum effect of several 
ions. Considering the AFC bandwidth T^pc — 60 MHz and the homogeneous linewidth of 
Nd 3+ ions Th = 63 kHz which have been measured under similar experimental conditions 



29| | . we can estimate that at least TAFc/^h ~ 10 3 different ions have to contribute to the 
AFC absorption and the envelope evolutions. The elementary AFC composed of about 10 3 
ions just likes a diffraction grating in the frequency domain, which governs the evolution 
of the atomic coherence and leads to a temporal delayed echo emission. One cannot fur- 
ther divide the elementary AFC grating into single ions to mimic the described evolution 
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FIG. 2: The experimental setup for violation of LGI with macroscopic crystals. The AOM1 
(AOM2) produces V (H) polarized pump light for the AFC preparation. The AOM3 produces 
a weak probe light for the storage sequence. The probe light's polarizations are controlled by 
the polarization beam splitters (PBS), the half- wave plate (HWP) and the phase plate (9). The 
samples are placed in a cryostat with a temperature of 1.5 K and a magnetic field of 0.3 T. Beyond 
the sample, a 45° HWP corrects polarization rotation caused by the sample. The probe light's 
polarizations are then analyzed with HWP and a Wollaston prism. The two MCs are used to 
protect the single-photon detector (SPD) from classical pump pulses. 

in our current experiment. The proposed test therefore goes beyond the demonstration of 
meso/macroscopic quantum coherence in collective phenomena. 
The crystals' state now can be represented by, 

(t) = 1/2(1 + e^'^^D) + 1/2(1 - e i{2 * 5t+ipo) )\A) (3) 

To probe the crystals' state, we need to collectively transfer the delocalized atomic state into 
a photon state. Thanks to the periodic structure of AFC, a strong rephasing echo occurs 
after a time t s = 1/A. The photon is re-emitted in the forward direction as a result of a 
collective interference between all of the atoms that are in phase. The crystal's state \D) 
and | A) corresponds to readout photon polarization of H + V and H — V, respectively. 

Fig. 2 shows the experimental setup for investigating the evolution of the macroscopic 
crystals. The laser source is a cw Tiisapphire laser (M Squared, Solstis). The H and 
V polarized pump light is independently generated with two 260 MHz AOM (Brimrose, 
RF bandwidth: 120 MHz) in double-pass configurations and combined with a polarization 
beam splitter (PBS). To protect the single-photon detectors (SPD) during the preparation 
procedure, two phase-locked mechanical choppers (MC) are placed in the pumping optical 
path and before the SPD, respectively. The H and V polarized pump light are mostly 
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absorbed by the first crystal and the second crystal, respectively. Carefully adjusting the H 
and V polarized pump power can optimize the storage efficiency and balance the efficiency 
of the retrieved H and V components. The photons to be stored are generated by another 
260 MHz AOM in double-pass configuration. The probe light are decreased to about 0.5 
photons per pulse by the neutral density filters (NDF). 

The pump light overlap with the probe light at the sample with a non co-linear config- 
uration. The probe light focuses to a diameter of about 100 [im, while the two pump light 
is independently collimated to produces a much larger diameter on the sample. The sample 
is placed in a cryostat (Oxford Instruments, SpectromagPT) at a temperature of 1.5 K and 
with a superconductor magnetic field of 0.3 T in the horizontal direction. The two parallel 
Nd 3+ :YV04 crystals' c-axes are placed in the horizontal direction. The probe light's polar- 
ization is controlled with a HWP and a phase plate (6). Because each crystal only strongly 
absorbs H polarized light, the H polarized input photon is stored in the first crystal, and 
the V polarized light is stored in the second crystal after polarization rotation by the HWP. 
To rotate the polarization back to the input state, another HWP is placed at 45° outside of 
the cryostat. The HWP and Wollaston Prism together choose the polarization of the single 
photon detections. The signal from the SPD is sent to the time-correlated single photon 
counting system (Ortec). 

The AFC preparation pulse sequence are shown in Fig. 1(b), in which the pump light's 
frequency is swept over 60 MHz in 240 fis and each frequency step has been assigned a 
specific amplitude to give an comb structure. The AOM1 and AOM2 are programmed to 
give two AFC structures with the same periodicity A and the frequency detuning of S. An 
experimental AFC structure with a periodicity of 20 MHz and 5 of 5 MHz is given in the 
Supplemental Information. The readout efficiency at different evolution time is shown in 
Fig. 1(c). We note that the AFC created directly in the frequency domain shows a better 
performance of readout efficiency and requires less laser power consumption. 

The pump and probe timing are controlled by two arbitrary function generators (AFG, 
Tektronix, AFG3252). To avoid the fluorescence noise caused by the classical pumping 
light, the probe cycle begins after waiting for time T w = 1.2 ms after the preparation cycle 
is completed. During the probe cycle, 1600 weak pulse with temporal width of about 25 ns 
are sent into the sample with a frequency of 1 MHz. The complete pump and probe cycles 
are repeated at a frequency of 40 Hz. A digital delay generator (SRS, DG645) provides the 
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FIG. 3: (a). Experimental test on the validity of stationary assumption. With an AFC prepared 
with a 5 of 5 MHz, the conditional probabilities QDA(t,t + 100 ns) (blue triangles), QDD(t,t + 33 
ns) (black squares), QAA(t,t + 67 ns) (red dos) and QAA(t,t + 100 ns) (green stars) show little 
dependence on the evolution time t. The Solid lines are the theoretical predictions, (b). With an 
AFC prepared with a 5 of 5 MHz and the initial state of D, the possibility to find crystals' state 
in D (black squares) and ^4(red dots) in different times, (c). The envelope evolution of K- (black 
squares) and K + (red dots). The Solid lines are the ideal quantum mechanical predictions and the 
blue dashed line is the classical bound. 

trigger signal for coincidence logic of single photon counting. 

We first present an experimental test on the validity of stationary assumption in our 
system. By controlling the probe photon's polarization, we can change the initial state of 
the crystals and let the crystals' state be D (A) at different readout time t (see methods). 
Fig. 3(a) shows the conditional probability Qij(t,t + r) to find a system in state j at time 
t + r, if it is in state i at time t. It can be seen that Qij(t,t + r) only depend on time 
difference r and show little dependence on the evolution time t within experimental errors, 
thus the stationary assumption holds in our system. 

We now analyze the LGI under the classical realistic description, where the two crystals 
can only be in one of the two states D and A. The input state is initially in p — D, 



S 



after evolution time t, the state becomes statistical mixtures pt = (1 — a)D + aA, where 
a represents the influence of the evolution with time t. We can get K(t) = 1 — 2a. With 
further identical evolution time t with the same process, the final state evolves to p2t — 
(a 2 + (1 - a) 2 )D + 2a(l - a) A. Therefore, Q(2t) = Aa 2 - 4a + 1. It is easy to verify that 
K(2t) - 2K(t) = Aa 2 > -1 and K(2t) + 2K(t) = 4(a - l) 2 > -1. 

Next, we analyzed the experiment from the viewpoint of quantum mechanics. Using 
Eq.(3) as the specification of our quantum state, we can easily obtain the correlation 
function K(t) = cos(27r5t). With the same analysis, K(2t) = cos(47r5t). The func- 
tions K± in Eqs.f TTP]) can be calculated as K_ = cos(47r<5t) — 2cos(27r<5t) and K + = 
cos(47r<5t) + 2 cos(27T<5t). It can be seen that K_ reaches its minimum value of -1.5 when 
2n5t = 7r/3 and K + reaches its minimum value of -1.5 when 27c5t = 2tt/3, which yields to 
maximal violation of the inequalities Eq.(l) and Eq.(2), respectively. 

Fig. 3(b) shows the corresponding evolution of the crystals' state with an initial state 
specified by D. With the AFC prepared with 5 of 5 MHz, the possibility to find the state 
D (A) oscillates with a periodicity of 200 ns. Based on the measurement results in Fig. 
3(b), the evolution of K_ and K + are shown in Fig. 3(c). The solid lines are theoretical 
predictions based on Eq. 3. K + reaches —1.44 ± 0.06 at 67 ns, which violates the classical 
limit of -1 by about 7 standard deviations. The value of K_ at 33 ns has been measured to 
be —1.34 ± 0.07, which deviates from the theoretical prediction. This is mostly caused by 
slow response of the AOM, which makes the imperfect separation of the transmitted photons 
and the readout photons for short storage time. 

The user-programmable AFC polarization detuning 5 determines the phase evolution 
speed of the crystals' state. The violation of LGI should hold for a different 5, i.e., a 
different evolution speed. We further show the envelope evolution of K_ and K + with a 5 
of 2.5 MHz in Fig. 4. The corresponding crystals' state evolution is shown in Fig. 4(a), 
it oscillates with a periodicity of 400 ns. K+ (K_) reaches -1.45 ± 0.06 (-1.43 ± 0.07) 
at 67 ns (125 ns), which violates the classical limit of -1 by about 7 standard deviations. 
We attribute the deviations from the theoretical predictions to the imperfect spatial mode 
matching of the retrieved H and V components. 

The physical system undergoing evolution can be further enlarged by choosing dopants 
with narrower homogenous linewidth. The figure of merit TAFc/^h can be larger than 10 7 
since a homogeneous linewidth of 50 Hz have been obtained in Er 3+ :YSO crystals [32J. In 
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(a) 




Time (ns) 

FIG. 4: (a). With an AFC prepared with a S of 2.5 MHz and the initial state of D, the possibility 
to find crystals' state in D (black squares) and A(red dots) in different times, (b). The enve- 
lope evolution of K_ (black squares) and K + (red dots). The Solid lines are the ideal quantum 
mechanical predictions and the blue dashed line is the classical bound. 

that case more individual ions will contribute to the elementary AFC gratings. The distance 
between the crystals is not limited to 2.4 mm as obtained in our experiment. Referring to 
another experimental realization of quantum memory for photonic polarization with similar 
setup [ijj], in which the optical distance between two Nd 3+ :YSO crystals is about 1 m. 
The high fidelity performance of that quantum memory demonstrates that increasing the 
distance between two crystals generally do not affect the coherence between them. Thus the 
two crystals can in principle be brought much further apart in our work. 

In summary, we experimentally observed the coherent superposition of delocalized states 
of macroscopic systems and controlled their evolution with time. These results demonstrate 
that the dynamics of a collective excitation which is distributed in two macroscopic crys- 
tals cannot described by a macrorealistic theory. Contrary to the situation encountered is 
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collective effects yielding to standard macroscopic quantum coherence, the phase relation 
between the considered states is governed by a macroscopic number of atoms and therefore 
the observed interference effects are macroscopic in origin although the interfering states 
themselves are of low Leggett-disconnectivity. Since the stationary assumption was exper- 
imentally verified in our system, the violation of this form of LGI represents a strict test 
against macroscopic realistic models of the described set up. These results confirm the per- 
sistence of quantum effects in systems of increasing complexity which can even be seen with 
naked human eyes. 



Methods 



Polarization dependent AFC preparation. The AFC preparation pulse sequence are 
shown in Fig. 1(b), it requires simultaneously sweeping the laser frequency and modulating 



its intensity. This pumping technique has been used previously in c 



. In our experiment 



the pump light's frequency is swept over 60 MHz in 240 [is and each frequency step has been 
assigned a specific amplitude to give an comb structure. This procedure ensures a spectral 
tailoring resolution of about 500 KHz. The pump sequences are continuously repeated in 
10 ms to achieve an optimal AFC structure. The duty ratio in each cycle is carefully 
adjusted to take into account the power broadening effect. The AOMs are all controlled by 
programmable RF drivers. The AOM1 and AOM2 are programmed to give a polarization 
dependent AFC structure with a user controlled frequency detuning of 5. Refer to the 
experimental AFC structure given in the Supplemental Information, we conclude that the 
more than twofold efficiency promotion compared to our previous experiment 11] is mainly 
caused by the unaffected peak absorption of AFC teeth. 

Experimental test on the validity of stationary assumption. The input photon's 
polarization is rotated to H + V (H — V) by the HWP. To detect the crystals' state in D (A) 
at different readout time t, the probe photon state should be set as H + e l,po V (H — e m> V) 
with y? = —2ii5t by carefully adjusting the phase plate 9. This photon state is directly 
mapped to the crystals' initial state. The subsequent system evolution is recorded as the 
function of the time difference r. The conditional probabilities Qij(t,t + r) only depends 
on the time difference r and have little difference for different time t, as shown in Fig. 3(a). 
These results demonstrate that the stationary assumption is valid in our system. 
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Supplementary Information 
Polarization dependent AFC structure 

In our experiment, an optimized efficiency of about 21% at 50 ns storage time is obtained. We 
measure the AFC structure by sending long probe pulse and measuring their transmission. 
An example AFC, with a periodicity of 20 MHz and 5 of 5 MHz, is shown in the Fig. SI. 
Compared with our previous paper, the efficiency promotion is mainly caused by the large 
available optical depth (~ 6.2). The peak absorption is preserved because the inhomogeneous 
absorption is directly tailored in frequency domain. The experimental efficiency fits well the 
theoretical estimation. For longer storage time, the readout efficiency drops because of a rising 
background absorption. 
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FIG. 5: An AFC prepared with a periodicity of 20 MHz and 5 of 5 MHz. The red circles are 
measured with H polarized probe light and the black squares are measured with V polarized probe 
light. 
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